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Abstract 



The method of analytic perturbation theory, which avoids the problem of 
ghost-pole type singularities and gives a self-consistent description of both 
spacelike and timelike regions, is applied to describe the "light" Adler function 
corresponding to the non-strange vector channel of the inclusive decay of the 
r lepton. The role of threshold effects is investigated. The behavior of the 
quark-antiquark system near threshold is described by using a new relativistic 
resummation factor. It is shown that the method proposed leads to good 
agreement with the "experimental" Adler function down to the lowest energy 
scale. 
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I. INTRODUCTION 



In studying the relationship between theoretical predictions and experimental data, it 
is important to connect measured quantities with the most elementary theoretical objects 
to check direct consequences of the theory without essential use of model assumptions. 
Some single-argument functions which are directly connected with experimentally measured 
quantities can play the role of these objects. A theoretical description of inclusive processes 
can be made in terms of functions of this sort. Among them is the Adler function |I| 
which can be extracted from the experimental data for the process of e~^e~ annihilation into 
hadrons and the inclusive decay of the r lepton. The mass of the r lepton, M^- = 1.777 GeV , 
is large enough in order to produce decays with a hadronic mode. At the same time, in the 
context of QCD, the mass is sufficiently small to allow one to investigate effects lying beyond 
the framework of the perturbative approach. At present, there is rich experimental material 
obtained from hadronic decays of the r lepton. The first theoretical analysis of hadronic 
decays of a heavy lepton was performed in 1971 well before the experimental discovery of 
the r-lepton in 1975. Since then, the properties of the r have been studied very intensively. 

The ratio of hadronic to leptonic widths for the inclusive decay of the r-lepton, 

^ r(''"^ hadrons Ur) 

is the most precise one for extracting of the values of the fundamental QCD parameters at 
a low energy scale 0. The initial theoretical expression for i?^ contains an integral over 
timelike momentum 

2 ds f s V f s \ 

which extends down to small s and cannot be directly calculated in the framework of 
standard perturbation theory (PT). Indeed, the hadronic correlation function n(s) is 
parametrized by the perturbative running coupling which has unphysical singularities and, 
therefore, is ill-defined in the region of small momenta. To avoid this problem, one usually 
applies the following procedure. The initial integral (Q) is rewritten by using the Cauchy 
theorem in the form of a contour integral in the complex plane with the contour running 
around a circle with radius 

where D{z) = —z dll{z) / dz is the Adler function. This trick allows one, in principle, to avoid 
the problem of a direct calculation of the Rr ratio by Eq. ([T|). However, in order to perform 
this transformation self-consistently, it is necessary to maintain correct analytic properties 
of the hadronic correlation function, which are violated in the framework of standard PT. 
The analytic approach to QCD 0, the so-called analytic perturbation theory (APT) 
maintains needed analytic properties and allows one to give meaning to the initial expression. 
The APT description can be equivalently phrased either on the basis of the expression (|1|) or 
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on the contour representation (0) |^. The information obtained in r measurements allows 
one to construct various "experimental" curves. In particular, in paper a quantity i?^(so) 
with a variable "mass" Sq < M^, has been considered. This quantity, defined for timelike 
momenta, is similar to the "smeared" functions constructed, for example, according to the 
Poggio-Quinn- Weinberg method ||T0| . The r-decay data allow us also to determine quantities 
in the Euclidean region, including the Adler function. This function can be extracted from 
the r data down to the lowest energy scale |TT] . 



In this paper we study the Adler D-function corresponding to r decay mediated by the 
non-strange vector current via W~ — > du. We use the APT method which does not encounter 
the problem of unphysical singularities of the running coupling and gives a self-consistent 
description of both the timelike, Eq. ([l|), and the spacelike, Eq. (H), regions. 

The region of integration in Eq. (|1|) includes the vicinity of the quark-antiquark threshold. 
The perturbative expansion breaks down in this neighborhood due to singularities at s = 
(mg -|- rrig)'^ |]T^,|TD|]. Thus any finite order of the perturbative expansion is unreliable near 
quark thresholds and, therefore, all singular terms of the (a^/f)" type, where v is the 
relative velocity of the quarks, have to be summed. Note that this problem cannot be 
avoided by using the contour representation (||) instead of Eq. (|l|), because these expressions 
should be equivalent to each other in the framework of a consistent method. For heavy 
quark systems one usually uses the nonrelativistic resummation factor obtained by using 
the Schrodinger equation with the Coulomb potential, which is known as the Sommerfeld- 
Sakharov factor But for a systematic description of the threshold region in the 

system of light quarks it is necessary to apply a relativistic approach. Here, we take into 



account threshold effects by using a new relativistic resummation factor proposed in Ref. []T5 
which was obtained for a QCD-like potential. 

II. ANALYTIC APPROACH TO r DECAY 

We start our consideration with a three-loop PT and APT analysis, neglecting, in the 
beginning, quark masses. It is convenient to separate the QCD contribution by representing 
the Rr ratio in the form 

Rr = K{l + SQCB), (3) 

where -R° corresponds to the parton level description and 5qcd is the QCD correction. We 
introduce QCD contributions to the imaginary part of the hadronic correlator, r(s), and to 
the corresponding Adler function, d{z) as follows: 

n{s) = Im n(s + ie) /t^RI ocl + r, D{z)(xl + d{z) . 

Then, one can write 5qcd as an integral over timelike momentum (Minkowskian region) 

or as a contour integral in the complex plane (Euclidean region) 
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We now consider 5qcd in the framework the PT and APT methods. 

A. Perturbation theory 

The PT description is based on the contour representation and can be developed in the 
following two ways. In Braaten's (Br) method the quantity (^ is represented in the form 
of truncated power series with the expansion parameter a,- = as{M'^)/TT. In this case the 
three-loop representation for 5qcd is 

CcD = «r +rial + r2 , (6) 

where the coefficients ri and r2 in the MS scheme with three active flavors are ri = 5.2023 
and r2 = 26.366 §. 

The method proposed by Le Diberder and Pich (LP) uses the PT expansion of the 
d-f unction 

d{z) = a{z) + di a'^{z) + d2 a'^^z) , (7) 



where in the MS-scheme di = 1.640 and d2 = 6.371 |16] for three active quarks. The three- 
loop PT running coupling, a{z), is commonly written in the form of an expansion in inverse 
powers of L H . In the MS scheme it is 



4 1 fthiL 1 



ln(-z/A^) , (8) 



where /3o = 11 - 2n//3, /3i = 102 - 38n//3 and (3^^ = 2857/2 - 5033n//18 + 325^^/54 are 
the first three /3-function coefficients. The substitution of Eq. (0) into Eq. (||) leads to the 
following non-power representation 

^LP^ = A^^) (a) + di (a) + d2 ^(2) (a) (9) 

with 

A(-){a) = — i ^(i-'Y(i+']a^(,). (10) 

As noted above, transformation to the contour representation (^ requires the existence 
of certain analytic properties of the correlator: namely, it must be an analytic function in 
the complex ^-plane with a cut along the positive real axis. The correlator parametrized, 
as usual, by the PT running coupling does not have this virtue. Moreover, the conven- 
tional renormalization group method determines the running coupling in the spacelike re- 
gion, whereas the initial expression (||) contains an integration over timelike momentum, 
and there is the question of how to parametrize a quantity defined for timelike momentum 
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transfers |T^. To perform this procedure self-consistently, it is important to maintain correct 
analytic properties of the hadronic correlator [18-pO[|. Because of this failure of analyticity, 
Eqs. d^) and (|^) are not equivalent in the framework of PT and, if one remains within PT, 
it is difficult to estimate the errors introduced by this transformation. However, using the 
APT method, it is possible to resolve these problems.^ 



B. Analytic perturbation theory 

In the framework of the analytic approach,^ the functions d{z) and r{s) are expressed in 
terms of the effective spectral function p(cr) 

diz) = - r^pia), r{s) = - r^pia). (11) 
n Jo c — z n Js cr 

The APT spectral function is defined as the imaginary part of the perturbative approxima- 
tion to dpt on the physical cut 

pier) = eo{a) + diQiia) + d2e2icr) , (12) 

where 

gn{cr) = lm[a;+\a + te)]. (13) 

The function ^0(0") in Eq- dH) defines the analytic spacelike, aa.n{z), and timelike, aan(s), 
running couplings as follows 



Oar, Z 



°° da , , „ , , 1 /-^ rfa , , 



TT Jo a — Z TT Js a 



As has been argued from general principles, the behavior of these couplings cannot be the 
same [^. It should be stressed that, unlike the PT running coupling, the analytic running 
coupling has no unphysical singularities (the ghost pole and branch points) and, therefore, 
possesses the correct analytic properties, arising from Kallen-Lehmann analyticity refiecting 
the general principles of the theory. For example, the one-loop APT result is PJT9[] 



4 A- . 4 



1 1 ln(s/A 



2^ 



arctan 



2 TT TT 



(15) 



where (z) = 4/ [Po\n{-z/A^)] . 

The analytic running couplings (the exact two-loop and the three-loop after an approx- 
imation) can be written explicitly in the term of the Lambert function pSlES. However, 



^The nonperturbative a-expansion technique in QCD [21| also leads to a well-defined procedure 
of analytic continuation |18|. 



^To distinguish APT and PT cases, we will use subscripts "an" and "pt" . 
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in the framework of the APT approach there is a httle sensitivity to the approximation in 
solving the renormahzation group equation for the running couphng In the following, 

we use an explicit form for the analytic running coupling in the timelike region, d^nis), that 
is derived by using the formula (H) in Eq. (p!2[) . The corresponding expression reads 



s) + Ai(s) + Aafs) 



where a^^^{s) is given in Eq. (|T^) and 



Ai 
A2 



A{s] 



4vr/3f 1 J 

f3o 



\nB + l-A 



ln(s/A2 



TT 



PI 



Train (s) 



B(s) = 



ln(s/A2 



+ 



TT 



1 - 



7r2 + ln'(s/A2). 



AlnB 



(16) 



(17) 



For the QCD correction r(s) one can write down E5 



?^an(s) = aan(s) + Ari(s) + Ar2(s) 
where aan(s) is given by Eq. ([T6| ) and 
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1 



Pi TT 



Jo B^ 



PIB^ 



(21nE + l)MVA!)+^_^i!W 



TT 



Aro = ci2 



Jo B^ 



ln(s/A2) . 



(18) 



(19) 



Using Eq. (|D or equivalently Eq. (|^), we obtain the QCD correction to the Rr-TsXio in 
terms of p((j) as follows 



°° c/cr 



TT JO 



cr 



M2 



+ 



M2 



p(a). 



(20) 



It is obvious that the first term of this expression is ran evaluated at the r mass and that 
5an is not representable as a series expansion in the running coupling. 

The difference between the PT (LP) and APT contributions to Rr can be transparently 
shown by the one-loop relation: 



(21) 



The additional term, which is "invisible" in the perturbative expansion, turns out to be 
important numerically |l7|j26[]. Due to the negative sign of this term, the QCD scale is larger in 
this method, Aan > Apt, at the same value of the QCD correction: 5an(Aan) = ^pt(Apt) = 5"^^^ . 
It should be noted that due to the difference of shapes of the APT and PT running couplings 
(see, for example, Ref. [^), their values at the r scale do not differ very much . 
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The APT analysis of the r decay in the three-loop level has been performed in Ref. p8 . 



This investigation together with other results (see, for example, Refs. |^,29,30|) allows us 
to formulate the following features of the APT method: (i) this approach maintains the 
correct analytic properties and leads to a self-consistent procedure of analytic continuation 
from the spacelike to the timelike region; (ii) it has much improved convergence properties 
and turns out to be stable with respect to higher-loop corrections; (iii) renormalization 
scheme dependence of the results obtained within this method is reduced dramatically. 

C. The vector channel in r decay 

Experimentally the i?T--ratio can be separated into three parts 

Rj- = Rry + Rt,a + Rt,s ■ (22) 

The terms Rt,v and Rt,a are contributions coming from the non-strange hadronic decays as- 
sociated with vector (V") and axial- vector [A) quark currents respectively, and R^^s contains 
strange decays (5). 

Within the perturbative approximation with massless quarks the vector and axial- vector 
contributions to Rr coincide with each other 

3 

Rry = Rt,a = -\Vud\ (1 + ^qcd) , (23) 
where \Vud\ denotes the Cabibbo-Kobayashi-Maskawa matrix element. However, the exper- 



imental measurements |31,32| show that these components are not equal to each other. The 
corresponding difference is associated with non-perturbative QCD effects which are usually 
described in the form of power corrections. The experimental data for the isovector spectral 



function of the ALEPH Collaboration have been used in Ref. |TT]| to extract the Adler 
Dy-function which we show as the dashed line in Fig. |^. The experimental D-function turns 
out to be a smooth function without any trace of resonance structure. The D-function 



obtained in Ref. |^3[ from the data for electron-positron annihilation into hadrons also has 
a similar property. One can expect that the Adler D-function defined in the Euclidean 
region reflects more adequately the quark-hadron duality than do quantities determined in 
the Minkowskian region,^ and, therefore, is more suitable for relating theoretical predictions 
with experimental data. It is obvious that the PT approximation cannot be applied at the 
low energy scale. Besides, any finite order of the operator product expansion fails to describe 
the infrared tail of the D-function. APT's good behavior in the infrared region, in principle, 
allows us to consider the D-function down to a low energy scale. 



^The Minkowskian and Euclidean characteristics which are associated with the process of electron- 



positron annihilation into hadrons have been considered in Ref. [34|. 



7 



III. THE VECTOR CHANNEL L»-FUNCTION 



In this section we describe the Adler function corresponding to inclusive r decay. By 
using a dispersion relation for the hadronic correlator n(g^), one can represent the Adler 
function as follows 

We will derive TZ{s) by applying the APT approach while incorporating threshold re- 
summation. 



A. Threshold effects 



A description of quark- ant iquark systems near threshold requires us to take into ac- 
count the resummation factor which summarizes the threshold singularities of the pertur- 
bative series of the [as/v)^ type. In a nonrelativistic approximation, this is the well known 
Sommerfeld-Sakharov factor [0,0 . For a systematic relativistic analysis of quark- ant iquark 



systems, it is essential from the very beginning to have a relativistic generalization of this 
factor. Moreover, it is important to take into account the difference between the Coulomb 
potential in the case of QED and the quark- ant iquark potential in the case of QCD. This 
QCD relativistic factor has been proposed in Ref. to have the form 



Six) 



1-exp l-Xix)] 



x{x) 



Att as 
3 sinh X 



(25) 



where x is the rapidity which related to s by 2m cosh x = V^- The relativistic resummation 
factor ( p5D reproduces both the expected nonrelativistic and ultrarelativistic limits and 
corresponds to a QCD-like quark- ant iquark potential. 

A convenient way to incorporate quark mass effects is to use an approximate expres- 
TO|JT^ which here can be written as 



sion 



n{s) = T{v) [1 + g{v)r{s)] Q{s - 4m^ 



where 



Tiv) 



4tt 



vr 
2^ 



3 + v 
4 
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Am? 



(26) 



(27) 



We introduce effective quark masses, which accumulate some non-perturbative contribu- 
tions and turn out to be close to the constituent masses. In the description of the non-strange 
vector component of the Z)-function we neglect the difference of the quark mass values and 
set m„ = rrid = m. 

The threshold resummation factor (^) leads to the following modification of the expres- 
sion (^) 



TZvis) = T{v) 



S{x)-\x{x)+g{v)r{s) 



6(s - 4m^ 



(28) 
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which one can use to calculate the vector component of the Rr ratio 



Rr.v = 3S,wlK/ I"' ^ (l - -ij) ' (l + -|) . (29) 



where Sew denotes the electroweak factor ||35| and the light vector /^-function is defined by 
Eq. (HD with n^TZv 



B. Results 

We derive the Z^vz-function and the value of Rry by using Eq. (|28|) with the expression 
(|I^) for r(s) and the timelike running coupling (|IBp in the S'-factor. In Fig. ^ we plot this 
Z)-function (solid curve) which was calculated by using the parameter A = 420 MeV and 
the value of the quark masses rriu = rrid = 250 MeV. Note that, practically, the same values 
of the light quark masses have been obtained in Refs. []36| , |37| . These values are close to the 
constituent quark masses and incorporate some nonperturbative effects. The shape of the 
infrared tail of the Z)-function is rather sensitive to the value of these masses. We obtain the 
value of Rry = 1.77 which agrees well with the experimental data presented by the ALEPH, 
^expt ^ ^ j^^ ^ Q Q^-^ OPAL, RlJ' = 1.764 ± 0.016 ig, collaborations. The 



value of the parameter A obtained here is slightly more than in the standard analysis. 
For example, the ALEPH Collaboration result is A = 370 ± IScxpt ± 38theor MeV. In this 
connection, note that a more accurate account of the so-called vr^-terms in the treatment of 
high energy experimental data for timelike processes leads to an increase in the value of as 
20| . It should be emphasized that the value of A turns out to be significantly smaller than 



the value extracted from the r data in Ref. p8|, which shows the importance of threshold 



effects. [Recall that in Fig. Q we also show the experimental curve (dashed line) and the 
curve corresponding to the perturbative result with power corrections (dotted line) which 
are taken from Ref. [111].] Note also that using the above values of quark masses, we find 
that the one-loop APT description gives a result (not shown in the Figure) which is close 
to the three-loop one and reproduces well the experimental behavior of the D-function. 



IV. CONCLUSION 

In this paper we have considered the Adler function corresponding to the non-strange 
vector channel data from r decay. This function, defined in the Euclidean region, is a 
smooth function and represents a convenient testing ground for theoretical methods. We 
have proposed the method of the "light" vector D-function description. The conventional 
method of approximating this function as a sum of perturbative terms and power correc- 
tions cannot describe the low energy scale region because both the logarithmic and power 
expansions diverge at small momenta. We have used the analytic approach to QCD which 
is not in conflict with the general principles of the theory and, in the infrared region, has a 
regular behavior. The new relativistic resummation factor, which corresponds to a QCD-like 
quark- ant iquark quasipotential, has been used to incorporate the threshold effects. We have 
shown that our approach allows us to describe well the experimental data for r-lepton decay 
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in terms of the D-function down to the lowest energy scale and for Rr in the non-strange 
vector channel. We have found that the influence of relativistic threshold resummation is 
important and leads to an significant reduction of the value of the QCD scale parameter A 
extracted from the r data. 
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FIG. 1. The light Adler function corresponding to the non-strange vector channel of r decay 
data. The solid curve is the APT result. The experimental curve (dashed line) corresponding to 
the ALEPH data and the perturbative result with power corrections (dotted line) are taken from 
Ref. Inl. 
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